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In numerical calculations of the primordial power spectrum of perturbations produced during 
inflation, for very little wavenumber ks, to implement the initial conditions required of the pertur- 
bation mode functions consistently, we must let the universe experience more e-foldings of inflation 
than that required of it to solve the horizon, flatness and other pre-inflation problems. However, if 
the number of e-foldings of inflation has an upper bound, then for perturbations at scales greater 
than some critical value, the initial conditions required of the mode functions cannot be imple- 
mented physically. Because at the inflation beginning point, these perturbation modes lie outside 
the horizon. We propose that such perturbations do not contribute to the Cosmological Microwave 
Background Anisotropy (CMBA). Under this proposition, the exceptional lowness of the observed 
little / muli-poles of CMBA is reproduced numerically. In Linde's 0^ model, the upper bound on 
the number of e-foldings of inflation is determined to be 65 approximately. 



I. INTRODUCTION 



this paper. 



Although calculations of the primordial power spec- 
trum of perturbations produced during inflation have 
been discussed by several authors 3, 6, 8], wc still think 
that there are problems worthy of notifying in this pa- 
per. For example, in the slow roll approximations, to get 
the nearly scale free primordial power spectrum of per- 
turbations, reference Q and Q employed the asymptotic 
expressions of the Hankel function as its argument goes 
to zero. This is difficult to understand because what we 
want to get is the value of the function when its argu- 
ment equals one. For details, please see our discussions 
in the second section. 

In the third section, we will discuss effects from the 
existence of an upper bound on the number of e-foldings 
of inflation on the implementation of the initial condi- 
tions required of the perturbation mode function. If the 
number of e-foldings of inflation has an upper bound, 
then comoving horizon of the universe will have an up- 
per bound as we trace back to the origin of the universe. 
Perturbations at scales greater than this bound will keep 
out of the horizon regardless how early we trace back in 
the history of the universe. By assuming that such per- 
turbations do not contribute to the angular power spec- 
trum of CMBA, we will find that the observed excep- 
tional lowness of the little / multi-poles of CMBA could 
be reproduced numerically. In Linde's (j)'^ model, the up- 
per bound on the number of e-foldings of inflation takes 
the value estimated by T. Banks and W. Fischler lOj. 

In our discussions, when specific inflation potentials 
are needed, we will take Linde's chaotic model 1^(0) = 
rn^cjp 12 as the illustrating examples. Except notified ex- 
plicitly, only density perturbations will be discussed in 
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II. CALCULATIONS OF THE PRIMORDIAL 
POWER SPECTRUM OF PERTURBATIONS 

Taking the simplest single fleld inflation as an exam- 
ple, to get the primordial power spectrum of perturba- 
tions from a given model and make decisive predictions 
about its effects on the large scale structure formation in 
the late time cosmology, we need to solve_the following 
coupled Friedmann-Mukhanov equations 
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As soon as eqs(^ are solved, spectrum of the primordial 
perturbations can be given as 
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Our conventions and formalisms in this paper are in 
agreement with those of 

In the conventional slow roll approximations, begin- 
ning from. 



_-H cj) 



z 



2a^H'^{l + e 



-e6 



2H 



(4) 



2 



taking e and 5 as constant, eqQ with the initial condition 
eq(|2l can be solved analytically, 
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In eqs® and Q, hI^\x) is the third class of Bessel func- 
tion or Hankel function of order v. H^{1) weakly depends 
on k through i/, \Hi,5{l)\ w 1.13. We note here that our 

expressions for P^{k), is a little different from that of the 
reference 3, 8] in the numerical factor. The authors there 
used the asymptotical expression of eq® as /cry ^ in 
calculating their primordial power spectrums, which is 
unnecessary and difficult to understand because at the 
point of horizon exiting, aH = k, kij 1 0. The same 
question occurs to the expressions of the gravitational 
wave perturbation power spectrum of reference Isl, 1^ . 

In practice, except in the case of power law inflation, in 
almost all the inflation models, the slow roll parameters 
Q vary with time. So, eqs©-® can only be taken 
as the first order approximations. To get more precise 
results, one convenient way is numerical calculation. In 
numerical calculations, eqslpQ) should be written in the 
following first order differential equations. 
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FIG. 1: Horizon axis, scale factor. Vertical axis, mode func- 
tion imaginary and real part, k = 0.3Mpc~^ . Uk os- 
cillates not only at the beginning of inflation but also at the 
end of it. The inflation driving potential is V((/f>) = m^(/!>^/2. 



In the above equations, we have taken scale factor a in- 
stead of time t or conformal time r as the independent 
variable. Reference ,6j and use Ina as the independent 
variable of numerical integration. In principle, using a or 
Ina is equivalent. But in practice, using a as the indepen- 
dent variable makes the stepsize control of the integration 
routines more convenient, so is more time saving. This 
is the first problem we think worthy of notifying in this 
note. The second problem is, for little(< 10~^Mpc~^) 
ks, we can integrate eas (llU|l both by the Fourth Or- 
der Runge-Kuta self-adaptive stepsize control integration 
method and the Blosch-Stoer method 0- But for large 
/cs, Runge-Kuta method becomes more and more time- 
costing so that it cannot adapt the variation of our in- 
tegration goal variables at all when k k, lAIpc~^. So, 
in our final CMBA multi-poles' calculation, we use the 
Blosch-Stoer method to calculate the primordial power 
spectrum of perturbations. 

For a given fc, we give our numerical results of the 
mode function Uk in FIG.Q From the figure we see that 
the mode function oscillates not only at the beginning of 
infiation, but also at the end of it. It oscillates at the 

end of inflation because Uk (77) (x ^ while </) has oscilla- 
tion behaviors during this time. However the quantity 
fe__|Mfcl!Z)| keeps stable after the first period decreasing 
during inflation, see FIG. ^ and 13 We do not know 
whether the oscillation behavior of Uk{r]) at the end of 
infiation will have any observation implications, we note 
it here only because it was ignored by peoples for long 
time in researches. 

In numeric calculations, we take the integration start- 
ing point as the conformal time zero point. Using the 
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oscillating properties (time translation invariant, up to a 
phase 0,01) of Mfc at the beginning of inflation, the initial 
conditions can be written as 



iv)-- 



1 



2k 
k 



-ik7] 



-ikr] 



(11) 



To solve the horizon, flatness and other pre-inflation 
problems, we require the universe expand Ne = 65 e- 
foldings during inflation. The initial values of the scalar 
field has been tuned so that 65 e-foldings of inflations 
can be obtained. Assuming the reheating temperature is 
about lO^^Gev and the universe expand adiabatically in 
the standard big bang process, then at the initial point 
of integration, the scale factor is 
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provided that the scale factor of today is set to 1 . How- 
ever, if we let k varies from [10'* Mpc]~^ to IMpc^^ and 
integrate eas ((Tn|) from a — 10~^^ to a = 10~^* to get the 
full power spectrum, what we find will be a large scale 
suppressed result, see FIG. |21 Inland the notes there. 

The observational result of jl|, if verified, favors a large 
scale suppressed power spectrum, but our suppressions 
here is due to an inconsistent initial conditions impos- 
ing. Because for very little fcs, at the integration starting 
point Gini = 10^^^, / {aH)~^- 0. So eqsdJ cannot 
be taken as the correct initial conditions. For example, 
if we want to integrate eqsJQ) in the simplest chaotic in- 
flation model, 1/(0) = m20V2, m = 5.61 x IQ-^Mpi for 
k = Hq = (4225Mpc)~*. To assure driving the universe 
expanding 10^* times during inflation, 0i„i = 16.2Mpi 
Then, at the integration starting point. 
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Obviously, k~^ / {aH)~^- ^ 0.60 ^ 0, i.e. this k mode 
of fluctuation is not well inside the horizon of that time. 
So, to use eas ((TT|) as the initial conditions, we have to 
integrate eqs([3l or equivalently eas H10|) beginning from 
a more little scale factor, say Umi = 10^^^ for example, 
so that those modes of fluctuations locate well inside the 
horizon of the integration starting point. Of course, the 
initial values of (p should be tuned correspondingly so that 
10'^" times of inflation can be obtained. Although works 
on the numerical calculations of the primordial power 
spectrum existed in several literatures 6, 7, 8], none of 
them pointed out this problem explicitly, so we think it 
is worth noting here. 
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during inflation, the problem discussed above would be 
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more serious. 



FIG. 2: Horizon axis: scale factor, normalized by setting 
cireh = 10~^*. Vertical axis: quantity -^^l^l- Different 
curves corresponds to different k. Short vertical lines in the 
figure cross the curves at aH — k point. Connecting centers 
of them yields the profile of the primordial power spectrum, 
see also the curves c and d of FIG.O Little figure: obtained 
by integrating eas llUII from a = 10~^® to a = 10~^*, taking 
egs im as the initial conditions. Large figure: the same as 
the little one, but the scale factor varies from 10~^* to 10~^* 
and the initial values of (j) is tuned correspondingly so that 
10'"' times of infiation can be obtained. The driving potential 
is 1/(0) = m^^VS- 



III. THE UPPER BOUND ON THE NUMBER 
OF E-FOLDINGS OF INFLATION 



However, if the number of e-foldings of inflation has 
an upper bound N^.b., then the comoving horizon of the 
universe will have an upper bound also. So as we trace 
back beyond this upper bound, the comoving horizon will 
decrease. Perturbations at scales larger than this upper 
bound will keep out of the horizon regardless how early 
we trace back to the origin of the universe. For such 
perturbations, the initial conditions eq® cannot be im- 
plemented physically. 

Then, what is the upper bound? Is it large enough so 
that, at least for perturbations at scales of today's Hub- 
ble horizon, the initial conditions eqs(|2Jl can be imple- 
mented numerically? In string theories, by T. Banks and 
W. Fischler's argument |[10|, if the present acceleration 
of the universe is due to an asymptotically de Sitter uni- 
verse with small cosmological constant, and the principle 
of Cosmological Complementary is valid, then the num- 
ber of e-foldings of inflation is bounded. They estimated 
that Nu.b. = 65. So the value already equals our as- 
sumed number of e-foldings of inflation required to solve 
the flatness, horizon and other pre-inflation problems in 
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FIG. 3: Primordial power spectrum of density perturbations 
in Linde's (f)^ model, c: no suppression at all. d: large scale 
power suppressed by the inconsistent initial conditions. 
1: by assuming that Pn{k) = for aU k < [7082Mpc]"^ 
by assuming that Pizik) = for aU k < [4225Mpc]"\ 
by assuming that Pn{k) = for aU k < [2113Mpc]"\ 



FIG. 4: Theoretical curves is computed by CMBfast. Ex- 
cept the primordial power spectrum, see FIG. Eland the notes 
there, other cosmic parameters are from WMAP's best fitting 
ACDM model. Scattered points are the observational results 
of WMAP2003 Q. 



eas (|13|l . Hence, if T. Banks and W. Fischler's estimation 
is the case, then at least for perturbations at scales of to- 
day's Hubble horizon, it will be impossible to implement 
the initial conditions physically. Hence the power 
spectrum of such perturbations could not be nearly scale 
free any more. It could be suppressed strongly comparing 
with that of the little scale perturbations. 

The second question is, if the number of e-foldings of 
inflation, consequently, the comoving horizon as we trace 
back to the origin of the universe has an upper bound, 
then, do perturbations at scales larger than this bound 
contribute to our observed CMBA? RecaUing ^ that, in 
calculating the multi-poles of CMBA, we need to do the 
following integration: 
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where P-jzik) is the primordial power spectrum of per- 
turbations produced during inflation and gj^{k) is the 
radiation transfer function. Mathematically, the integra- 
tion should be made from to cx3. But physically, if some 
mechanism exists so that perturbations at scales greater 
than a certain cutoff wavelength k^^ do not contribute 
to CMBA, then the lower bound of the integration l|14|l 
will be changed unavoidably. 

Obviously, if the number of e-foldings of inflation has 
an upper bound, then at the beginning of inflation, co- 
moving horizon of the universe will reach its upper bound 
[ai,b.Hi,i,]~^ as we trace back to the origin of the uni- 
verse. Perturbations beyond such a bound will keep out 



of the horizon regardless how early we trace back. We 
propose that such perturbations do not contribute to our 
observed CMBA. Our motivation is: from the beginning 
of inflation, these perturbations lie outside the horizon. 
They have no chances to experience a transition from 
quantum ones to classical ones before CMBA formed. 

Under this assumption, we write the primordial power 
spectrum of perturbations produced during inflation as 
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where A is a constant whose value can be determined 
by requiring that the corresponding primordial power 
spectrum gives approximately the same first peak CMBA 
multi-poles as observations, see FIG.|31and01 while fee is 
determined by the following relation: 



(16) 



where the subscript i.h. means that the relevant quantity 
take value when inflation begins. 
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In the above equations, the scale factor of today is set 
to one, Nu.b. denotes the upper bound on the number of 
e-foldings of inflation, Trek and Ttod are the temperature 
of the universe at the reheating point and today respec- 
tively, while Ainfi is the characteristic energy scale of a 
given inflation model. We supposed that the universe en- 
ter the standard adiabatic expansion immediately after 
inflation. 

Take T. Banks and W. Fischler's estimation as an in- 
put and eas ((T^ as a guide post, where by the logic of 
eqsdig N^,b. = 65, T,eh « lO^^Gev, Ttod ~ IQ-^ez;, 
Aj„/, = [m202/2]i/4 ^ 8 01 X IQ-^Mpi, we numerically 
calculated the CMBA multi-poles for three values of kc 

fcci,2,3 = [7082Mpc]-\ [A225Mpc]-\ [2113Mpc]-^ (18) 

where 7082Mpc could be considered as the largest sub- 
horizon perturbation scale under the assumption that 
Nu.b. — 65; 4225Afpc is the perturbation scale corre- 
sponding to our today's bubble horizon; while 2113Mpc, 
half of that. Our results are displayed in FIG. 13 and 
01 curves 1, 2, 3. From the figure we see that the little I 
multi-poles of CMBA is very sensitive to the cutoff of the 
primordial power spectrum. If kc takes a certain value 
between kc2 and kcs, then theoretical curves will tally 
with the observational results to a better degree. 

So, if the number of e-foldings of inflation has an upper 
bound, and by our proposition, perturbations outside the 
horizon at the inflation beginning point do not affect the 
large scale structure formation at late times, then the ex- 
ceptional lowness of the little I multi-poles of CMBA can 
be reproduced numerically. In Lindc's 0^ inflation model, 
which has characteristic energy A « 10^^ Mpi, the upper 
bound on the number of e-foldings of inflation is approxi- 
mately 65. In other inflation models, if the observational 
result are to be reproduced, by ea H17|) . more lower char- 
acteristic inflation energy will mean a more little upper 
bound on the number of e-foldings, and vice versa. 

When this paper is posted in the network, we are noted 
that in reference the authors studied the stochastic 



inflation mechanism and got similar results like ours. In 
some sense, the mechanisms there can be assimilated to 
another statement of a limited number of e-foldings of 
inflation. 

IV. SUMMARY 

In this paper, we first present four notes on calcula- 
tions of the primordial power spectrum of perturbations 
produced during inflation. 1st, to get the nearly scale 
free power spectrum of perturbation, it is unnecessary to 
use the asymptotic behavior of the Hankcl function as its 
argument goes to 0. 2nd, scale factor a is a better inde- 
pendent variable than Ina in the numerical integrations 
of the Friedmann-Mukhanov equations. 3rd, the pertur- 
bation mode function Uk oscillates both at the initial and 
the end of inflation times. 4th, numerically, to implement 
the initial conditions of the perturbation equations con- 
sistently, we must let the universe experience more num- 
ber of e-foldings inflation than that required of it to solve 
the horizon, flatness and other pre-inflation problems. 

We then explore effects from the existence of an up- 
per bound on the number of e-foldings of inflation on 
the implementation of the initial conditions required of 
the perturbation mode function. When such an upper 
bound exists, the comoving horizon of the universe will 
have an upper bound also as we trace back to the origin 
of the universe. So perturbations at scales beyond this 
bound will keep outside the horizon regardless how early 
we trace back in the history of the universe. After propos- 
ing that such perturbations do not affect the formation 
of CMBA, we find the observed exceptional lowness of 
the little I multi-poles of CMBA could be reproduced 
numerically. In Linde's 0^ model, the upper bound on 
the number of e-foldings of infiation is determined to be 
65 approximately. 
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